Platoon formation in a traffic model with variable anticipation 
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A cellular automaton model of traffic flow taking into account velocity anticipation is introduced. 
The strength of anticipation can be varied which allows to describe different driving schemes. We 
find phase separation into a free-flow regime and a so-called t;-platoon in an intermediate density 
regime. In a w-platoon all cars move with velocity v and have vanishing headway. The velocity 
v of a platoon only depends on the strength of anticipation. At high densities, a congested state 
characterized by the coexistence of a 0-platoon with several u-platoons is reached. The results are not 
only relevant for automated highway systems, but also help to elucidate the effects of anticipation 
that play an essential role in realistic traffic models. From a physics point of view the model is 
interesting because it exhibits phase separation with a condensed phase in which particles move 
coherently with finite velocity coexisting with either a non-condensed (free-flow) phase or another 
condensed phase that is non-moving. 

PACS numbers: 45.70.Vn, 02.50.Ey, 05.40.-a 



I. INTRODUCTION 

In the last years, the continuous increase of traffic de- 
mand has prompted authorities around the world to place 
more emphasis on improving the efficiency and capacity 
on the roadway systems. Ecological considerations, space 
and budgetary constraints have limited solving traffic 
congestion by upgrading and constructing new roadway 
systems. Advanced technologies for vehicular traffic have 
been developed as a mean to improve the management 
of existing system and thus to solve traffic congestion, 
environmental issues and improve traffic safety. How- 
ever, to achieve these aims, an accurate forecast of the 
impact of these technologies is critical before their final 
deployment. 

Testing these advanced technologies on real traffic is 
not always feasible. In contrast, computer simulations as 
means for evaluating control and management strategies 
in traffic systems have gained considerable importance 
because of the possibility of taking into account the dy- 
namical aspects of traffic (see for example, 0,0, Q) and 
assessing the performance of a given advanced technology 
in a short time. 

Cellular Automata (CA) models for traffic flow 0, Q 
have shown the ability to capture the basic phenomena 
in traffic flow |fj. Cellular automata are dynamic mod- 
els in which space, time and state variables are discrete. 
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Discrete space consists of a regular grid of cells, each one 
of which can be in one of a finite number k of possible 
states. All cells are updated in discrete time-steps. The 
new state of a cell is determined by the actual state of the 
cell itself and its neighbor cells. This local interaction al- 
lows us to capture micro-level dynamics and propagates 
them to macro- level behavior. The discrete nature of 
CA makes it possible to simulate large realistic traffic 
networks using a microscopic model faster than real time 
0- 0- 13 • Now, almost eleven years after the introduction 
of the first CA models, several theoretical studies and 
practical applications have improved the understanding 
of empirical traffic phenomena (see e.g. 0, 0, 0, ITl|). 
Moreover, CA models have proved to be a realistic de- 
scription of vehicular traffic in dense networks 0, • 

In this paper, we propose a single-lane probabilis- 
tic model based on the first CA model of Nagel and 
Schreckenberg Q (hereafter cited as the NaSch model) 
to describe the effects of several anticipation schemes in 
traffic flow. Anticipation in traffic means that drivers 
estimate their leader's velocities for future time-steps 
[H [H Q E E [H El. This can lead to an in- 
crease of the vehicular capacity and a decrease of the 
speed variance. However, incorporating different driv- 
ing strategies requires a safety distance with respect to 
the preceding vehicles. For that purpose we introduce a 
new parameter in the deceleration process, called antici- 
patory driving parameter, to estimate the velocity of the 
precedent vehicle. This estimation, plus the real spatial 
distance to the leading vehicle, determine a safe braking 
distance. By appropriately tuning this new parameter 
different traffic situations of non-automated, automated 
and mixed traffic can be considered. Furthermore, the 
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anticipatory driving parameter is relevant for Automated 
Highway Systems (AHS) [Hill. 

According to simulation results from our proposed 
model, the relations derived from the density vs. ve- 
locity and density vs. flow curves are in agreement with 
the fundamental diagrams that describe these relations 
in real non-automated traffic. In addition, simulation re- 
sults from our model in the case of high anticipation (like 
automation) describe one of the interesting phenomena in 
traffic flow, formation of platoons. We will show that, in 
contrast to models without anticipation, dense platoons 
can be formed where all cars move coherently with some 
finite velocity v > 0. The mechanism for platoon forma- 
tion is not only of great importance for AHS to increase 
highway capacity in a much safer way |2*lj , but also helps 
to understand the essential role of anticipation effects in 
realistic traffic models. By varying the anticipatory driv- 
ing parameter three different regimes, characterized by 
different slopes of the fundamental diagrams, can be ob- 
served. Apart from a free-flow and a congested phase, 
an additional regime where platoons of cars moving with 
the same velocity v < v max — 1 exist. 

The paper is organized as follows. In Section [n] we de- 
fine a modified NaSch model to consider different driving 
strategies. In this way, we introduce a new parameter 
to determine a velocity-dependent braking distance. In 
Section EU we present the results of our investigations. 
We show results for the fundamental diagram and differ- 
ent values of the anticipation parameter. A comparison 
with those for other models and real non-automated traf- 
fic is presented. Phase separation into a free-flow regime 
and so-called f-platoon is observed in a certain interme- 
diate density regime. For large densities, in the congested 
state phase separation into a dense jam (0-platoon) and 
w-platoons is observed. The flow structure determined 
by the existence of dense platoons with velocity v is cal- 
culated. Dependence of an optimal anticipation level on 
the density is found. Analytical results are in excellent 
agreement with results from computer simulations. In 
the concluding Section ll VI we summarize our results and 
discuss the relevance of our results for traffic models and 
real traffic. 



II. DEFINITION OF THE MODEL 

In this section we introduce the proposed model. It is 
defined on a one-dimensional lattice of L cells with pe- 
riodic boundary conditions, which corresponds to a ring 
topology with the number of vehicles preserved. Each cell 
is either empty, or it is occupied by just one vehicle trav- 
eling with a discrete velocity v at a given instant of time. 
All vehicles have a velocity that ranges from 0, . . . , v max . 
In addition, and for simplicity, only one type of vehicles 
is considered. The time-step (At) is taken to be one sec- 
ond, therefore transitions are from t — > t + 1. It can also 
be easily modified. 

Let Vi and Xi denote the current velocity and position, 



respectively, of the vehicle i, and v p and x p be the velocity 
and position, respectively, of the vehicle ahead (preceding 
vehicle) at a fixed time; di := x p — Xi — I denotes the 
distance (number of empty cells) in front of the vehicle 
in position Xi that sometimes is called headway. 

The dynamics of the model are defined by the following 
set of rules, that are applied to all N vehicles on the 
lattice each time-step: 

Rl: Acceleration 

If Vi < v max , the velocity of the car i is increased 
by one, i.e., 

vi -> mfn(Vi + l,v max ). 

R2: Randomization 

If Vi > 0, the velocity of car i is decreased randomly 
by one unit with probability R, i.e., 

Vi — > max(vj — 1,0) with probability R. 

R3: Deceleration 

If d^ < Vi, where 



(l-a)-v p + 



with a parameter < a < 1, the velocity of car i 
is reduced to df. [x] denotes the integer part of x, 
i.e. [x + i] corresponds to rounding x to the next 
integer value. 

The new velocity of the vehicle i is therefore 
Vi — > min(wj, df). 

R4: Vehicle movement 

Each car is moving forward according to its new 
velocity determined in steps 1-3, i.e., 

Xi ► Xi "T Vi . 

Rules Rl, R2 and R3 are designed to update velocity of 
vehicles; rule i?4 updates position. According to this, 
state updating is divided into two stages, first veloc- 
ity, second position. Note that this division follows the 
scheme in differential equation integration that first up- 
dates the time derivative and then the value of the state. 
It is important to mention that we are changing the order 
of the rules in comparison with NaSch model since R2 is 
applied before i?3. 

Rule Rl indicates that all the drivers would like to 
reach the maximum velocity when possible. Rule R2 
takes into account the different behavioral patterns of 
the individual drivers in which with no apparent reason 
a driver decreases its speed. These situations include, 
for example, cases of overreaction in braking or incidents 
along the highway that distract drivers, and random fluc- 
tuations. 
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Rule i?3 is the main modification to the original NaSch 
model Q . In this rule the distances between the ith and 
(i + \)th vehicles, and their corresponding velocities are 
considered. Knowledge of the preceding vehicle's velocity 
is incorporated through the anticipatory driving parame- 
ter a with range < a < 1. Notice that, by only varying 
the parameter a in the term df = di + [(1 — a)v p + 1/2], 
different anticipatory driving schemes that require differ- 
ent safe braking distance with respect to the preceding 
vehicle can be modelled. If a takes its maximum value 
(a = 1) the speed of the vehicle ahead is not considered 
in the deceleration process, i.e. anticipation is not consid- 
ered. On the contrary, when a = the speed of the vehi- 
cle ahead is considered without restrictions, i.e., without 
establishing a braking distance with respect to the prece- 
dent vehicle [2(j ■ This last case occurs with either a very 
aggressive driver or when vehicles can obtain information 
about the velocity of vehicles ahead 29] to allow small 
distances between vehicles (e.g. of the order 1 m). In- 
termediate values for a thus represent different braking 
spacing policies or degrees of automation in the vehicles 
or anticipatory driving schemes. Platooning schemes pH 
imply values of a closer to zero and demand additional 
req uirements to preserve safety, like coordinated braking 
19]. Independent vehicle driving with low level of antici- 
pation implies values of a closer to 1 in order to preserve 
safety levels: the larger a is the larger braking distance 
is. Note that in order to determine v n consistently for 
all vehicles in the case of periodic boundary conditions, 
ride i?3 must be iterated at most {v ma x — 1) times. In 
real situations, the drivers always estimate the velocity of 
preceding vehicle and according to this and their way to 
drive (relaxed or aggressive behavior) they choose a safe 
headway distance to drive. Variation of a allows also to 
model these aspects. 

Thus, the proposed model is able to represent different 
anticipatory driving schemes, and model the minimum 
braking distance required with only one parameter a, 
here referred to as anticipatory driving parameter. 

We emphasize that the CA model as presented here is a 
minimal model in the sense that all four steps BA-RA are 
necessary to reproduce the basic features of real traffic, 
however, additional rules may be needed to capture more 
complex situations 13]. 



III. SIMULATION RESULTS 

To simulate the CA model proposed in the previous 
section, the typical length of a cell is around 7.5 m. It is 
interpreted as the length of a vehicle plus the distance be- 
tween cars in a dense jam, but it can be suitably adjusted 
according to the problem under consideration. With this 
value of the cell size and a time-step of 1 s, v = 1 corre- 
sponds to moving from one cell to the downstream neigh- 
bor cell in one time-step, and translates to 27 km/h in 
real units. The maximum velocity is set to v ma x = 5, 
equivalent to 135 km/h. The total number of cells is as- 



sumed to be L — 10 4 , and the density p is defined as 
p = N/L, where N is the number of cars on the highway. 
Initially, N vehicles are distributed randomly on the lane 
around the loop with an initial speed taking a discrete 
random value between and v max . Since the system is 
closed, the density remains constant with time. 

Velocities are updated according to the velocity updat- 
ing rules Rl — R2 — R3 and then all of cars are moved 
forward in step i?4. Each run is simulated for T = 6L 
time-steps. In order to analyze results, the first half of 
the simulation is discarded to let transients die out and 
the system reach its steady state. For each simulation a 
value for parameter a is established by taking into ac- 
count the desired anticipation degree and thus, control- 
ling the safe braking distance among vehicles. In the 
following, the value of a is the same for all vehicles (ho- 
mogeneous drivers). 

A. Comparison with real non-automated traffic 

The fundamental diagram is one of the most important 
criteria to show that the model reproduces traffic flow 
behavior. This diagram characterizes the dependence of 
the vehicles flow on density. We obtain a fundamental 
diagram for the proposed model with R = 0.2 and a = 
0.75, see Fig. ^ This a value corresponds to cautious 
estimation of the preceding car's velocity. 

As we can see from Fig. ^ the curve of the proposed 
model (dashed line) is consistent with the characteristic 
curve of the measured fundamental diagram (solid line) 
taken from |23j |. The critical density and the maximum 
flux of our model are (p c ,q m ax) = (16%, 2417 cars/h), 
closer to the empirical curve values of (p c , q ma x) = (17%, 
2340 cars/h), in comparison with other existing models, 
see Note the quantitative agreement between simu- 
lated and empirical curve of the fundamental diagram in 
its decreasing part. 

Moreover, as we can see from Fig. Q] a decrease 
in the experimental curve slope is observed at (10%, 
1692 cars/h), corresponding to a reduced mean veloc- 
ity of vehicles near the critical density. Thus, we make 
a small modification to the proposed model based on an 
idea of |24| . We modified the deceleration rule such that 
a vehicle can only reach the maximum velocity if its brak- 
ing distance to the precedent vehicle is greater than nine 
cells (67.5 m, i.e, the density is lower than 0.1). The 
modified rule R2> of the model is then as follows: 

R3': If Vi = 5 (i.e. v.- L = v max ) and df < 9 cells then 
Vi — » mm(vi — 1, <if) 
else (like rule R3 of our original model) 
Vi — > min(ui,df). 

In Fig. |21 the results obtained from our modified model 
are shown and compared with the experimental curve 
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B. Modeling different anticipation schemes 




ocupancy (%) 

FIG. 1: Typical form of an empirical fundamental diagram 
taken from 23] (solid line) in comparison with simulations of 
the proposed model for R = 0.2 and a — 0.75 (dashed line). 



Determination of the impact of different driving strate- 
gies is important in order to propose automated traffic 
alternatives. Following that proposal, we decided to in- 
vestigate traffic flow behavior using our model. As men- 
tioned above, the parameter a represents the way in 
which different driving strategies adopt a braking dis- 
tance with respect to the preceding vehicles. Varying 
the parameter a, these strategies can be tuned. In Fig. El 
we show the fundamental diagram of the proposed model 
with a fixed value of R = 0.2 and different values of a. 
From this diagram, the impact of the driving strategies 
coded in a can be observed. Smaller values of a imply 
larger flows, that is, higher levels of anticipation. Here 
vehicles keep a less safe braking distance, leading to an 
increase in the vehicular capacity. This behavior is in 
agreement with, for example, platooning strategies that 
exploit the knowledge of the velocity of precedent vehi- 
cles and require a smaller distance among vehicles (near 
I m), so increasing the flow. 




.0 



occupancy[%] 



FIG. 2: Fundamental diagram for the modified model (R3'). 

This diagram is in agreement with the experimental curve in FIG - 3: Fundamental diagram for different values of the an- 
a quantitative way (R = 0.2 and a = 0.75). ticipation parameter a and R = 0.2. 



from real traffic data [23j. The curve resulting from the 
modified model fits the real curve in both its increasing 
part and its decreasing part quite well. Thus, the results 
from the modified model agree quantitatively with the 
experimental shape of the fundamental diagram. This 
is an indication that in real traffic the value of a may 
depend on velocity. 

In this subsection, we have shown that the given model 
can reproduce some common characteristics of the real 
manual (non-automated) traffic flow. However, with the 
new parameter it is possible to consider several antici- 
pation schemes for traffic flow: non-automated, mixed 
and automated traffic flow. This will be shown in the 
following for the original model (R3) defined in Sec. HTl 



It is important to notice that for values of a from 0.13 
to 0.50 a second positive slope corresponding to a mixed 
branch is observed in the fundamental diagram. It is in- 
teresting that the initial positive slope, corresponding to 
a free-flow region where there are no slow vehicles, is sim- 
ilar for all values of a. Here the vehicles travel at near 
maximum speed. For the second branch, on the other 
hand, the flow is increased with non-maximum veloc- 
ity, indicating a mixed region due to anticipation effects 
(Fig. BJ . In order to analyze the role of the anticipation, 
we show the average velocity as function of the density 
for the same parameter values as in Fig. |3J 

As we can see from Fig.0J higher levels of anticipation 
(smaller values of a) imply a larger density interval for 
the free-flow region. For values 0.13 < a < 0.50, after 
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0.0 0.2 0.4 0.6 0.8 1.0 

density 

FIG. 4: Relationship between mean velocity and density for 
R = 0.2 and different values of a. 



the free-flow region, traffic flow organizes in a so-called 
mixed region with a lower mean velocity. In this mixed 
region, in addition to free flowing vehicles, and vehicles 
moving in platoons where all cars have the same velocity 
and vanishing headway exist. The existence of this mixed 
region indicates that a suitable estimation of the velocity 
of precedent vehicle, coded in a, allows that more cars 
to fit on the road and the flow increases ever for values 
of large-density. 




FIG. 5: Space-time diagram showing the time evolution of a 
system simulated initially with a — 0.51. After switching to 
a — 0.5, the behavior changes dramatically. 

Fig-Elshows a space-time diagram for R = 0.2, p = 0.5 
that exemplifies the dramatic changes in the microscopic 
structure when changing the value of a |3dj |. The sim- 
ulation is started with a = 0.51 where the jamming re- 
gions travel backwards. However, after some time we 
switch to a = 0.5 and immediately observe a dramatic 



change in the slope of the congested regions. They are 
now traveling forward. This behavior has been observed 
before in anticipatory modeling 18]. Such structured 
flow observed in space-time diagrams increases the high- 
way capacity due to the space reduction among vehicles. 
In these simulations we have found that the branches of 
congested or jammed flow collapse to a single region as in 
the VDR model 25] . These results are analogous to those 
for slow-to-start models, because effectively the outflow 
from a jam is reduced compared to the maximal flow. 
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density 

FIG. 6: Standard deviation of the speed. 

On the other hand, it is also important to analyze the 
efficiency of traffic. In non-equilibrium situations, this is 
usually evaluated by analyzing the entropy of the system. 
However, for non-physical systems the analysis of the 
standard deviation of the main variable can give equiva- 
lent information j2|| . In the context of traffic flow, high 
standard deviation of speed means that, on average, a 
vehicle would experience frequent speed changes per trip 
through the system. In turn, the high speed variance 
could also increase the probability of traffic accidents. 
Therefore, the standard deviation of speed can be seen 
as an indicator of the efficiency in traffic flow. 

In order to analyze the efficiency of traffic with dif- 
ferent automation levels, we calculate the standard de- 
viation of speed (Fig. [fjj). For each value of a a maxi- 
mum that occurs shortly after the free-flow region can be 
clearly seen. In the free-flow region, the speed variance is 
negligible since there are no slow vehicles and fluctuations 
are extremely rare. Since the free-flow region increases 
as a decreases, it is seems reasonable to attempt traffic 
with the higher level of anticipation in the range of den- 
sity from to 0.5. This selection not only produces a 
state with higher flow, but also the lowest standard de- 
viation, so the efficiency is the highest. It can be clearly 
seen from Fig. that beyond the efficient density range, 
the level of anticipation coded in a should be switched 
based on the density regime: for p e (0.5,0.54], an effi- 
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cient performance is found with a — 0.13; however, for 
p € (0.54, 0.63] the highest efficiency is attained with 
a = 0.20. Summarizing, the behavior observed in Fig. 
indicates that the automation level should be determined 
depending on the density: higher densities require a lower 
automation level, that is, a safer braking distance. 

Besides, we can also see from Figure® the standard de- 
viation of speed resulting from assigning a random value 
of a between and 1 to each vehicle. We stress that 
for each car, this value of a is not changed during the 
time evolution. Therefore, we simulate, by means of ran- 
dom values of a, traffic with non-homogeneous drivers: 
aggressive, non-aggressive and relaxed drivers. We can 
see from Fig. [5] the region where the standard deviation 
is negligible is close to that resulting from homogeneous 
drivers with a = 0.2 (high-level anticipation). However, 
the variance of vehicles speeds is 60% larger than that 
corresponding to homogeneous drivers. 

In order to elucidate the effects of anticipation with 
non-homogeneous drivers, the fundamental diagram ob- 
tained for this random a case is compared with that 
for a = 0.2 in a homogeneous system (Fig. [7J|. Note 
that the second slope corresponding to a mixed region 
is missing in the inhomogeneous system. In this case, 
the variance in the level of anticipation considered by 
the drivers produces higher fluctuations of speeds, and 
the flow decreases rapidly. Therefore some anticipation 
driving schemes have a strong impact on the efficiency of 
the system. All of these findings require physical expla- 
nations |27j . 



herently with the same velocity v. In the following these 
will be denoted as w-platoons. The stationary state then 
shows phase separation into a free-flow region and a v- 
platoon. This is similar to the behavior observed in mod- 
els with slow-to-start rules where the system separates 
into free flow and a dense jam, i.e. a 0-platoon |25j. 

Fig. [S] shows the velocity distributions of the different 
branches. It can be clearly seen that only cars with veloc- 
ity v and free flowing cars (with velocity v max or v max — 1 
due to the randomization) exist. 



0.8 -i 
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FIG. 8: Velocity distributions for density p = 0.4 and different 
values of a. 
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FIG. 7: Comparison of the fundamental diagrams of homoge- 
neous systems with the random a case. 



C. Structure of the mixed states 

The behavior in the mixed states is determined by the 
existence of dense platoons in which vehicles move co- 



Since the headway d% of a car i inside a u-platoon is 
di = 0, its new velocity is determined by v[ — min(Uj, df) 



with df = [(1 - a) 



must be equal to v. 
condition: 



] . For a stable w-platoon, 



This gives the following stability 



v < (1 — a) • v + 



(1) 



Equation JTJ can be regarded as a condition for the an- 
ticipation parameter a. It implies that a w-platoon can 
only be stable for 



a v +i < a < a v 



where a v is defined by 



1 

2v' 



(2) 



(3) 



However, this condition is only necessary, not sufficient. 
The u-platoons that can be realized for a given a also 
depend on the randomization R. E.g. for R = 0.2, pla- 
toons with v = 0,1,2,3 occur, whereas for R = 0.4 
platoons with v = 3 can not be observed in the in- 
finite system although they might exist in small sys- 
tems. Simulations indicate that the slope of the mixed 
branch in the fundamental diagram has to be smaller 
than (1 — R)(vf — 1) where Vf is the average velocity 
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in free flow: Vf = v max — i?. This will be discussed in 
Sec. 1111 El in more detail. 

Another criterion for the stability of platoons can be 
derived from the condition that the inflow and outflow 
of the platoon have to be identical in the steady state. 
In the following we will derive estimates for these flows 
and in this way obtain the fundamental diagram in the 
mixed region. 

The outflow from a w-platoon is determined by the 
average time T w needed by the leading vehicle of the pla- 
toon to accelerate to velocity v + 1. Assuming that this 
car has a large headway, this time is determined by the 
randomization constant R through T w = j-rji- There- 
fore, in the free-flow region of the system, the average 
headway Axf is given by Axf — T w (yf — v) + 1. This 
consideration is very similar to the reasoning used in |25| . 

Assuming that the platoon consists of N v and the free- 
flow region of Nf vehicles, we have 



N = N v + N f , 



and 



L = N v + NfAxf. (4) 



Here N is the total number of cars. Furthermore it has 
been assumed that the transition region between the pla- 
toon, where all cars have headway di = 0, and the free- 
flow region, where the average headway is given by Axf, 
can be neglected. Eliminating Nf we find 



Nv_ 
L 



pAxj 



1 



Ax 



f 



1 



(5) 



We now can calculated the flow J = pv of the corre- 
sponding phase-separated state. The average velocity v 
in the presence of a v-platoon is given by 



N v v + NfV f 
N ' 



(0) 



A straightforward calculation using the results given 
above yields for the flow 



J v = {l-R) + {v-(l-R))p. 



(7) 



These results are in excellent agreement with the re- 
sults from computer simulations. Since 1 — R < 1 it 
implies that all slopes corresponding to mixed states in 
the fundamental diagram are positive, except those for 
0-platoons (a > 0.5 see Fig. |3J) which are responsible for 
the jammed branch with low flow and negative slope. 



D. Structure of the congested states 

For large densities all anticipatory curves (a < 0.5) col- 
lapse on one congested branch where the flow decreases 
with increasing density. Simulations indicate that the 
structure of the corresponding states depends on the pa- 
rameter regime. In the range PJl. where a v-platoon can 
exist, the congested branch is characterized by the co- 
existence of a compact jam (0-platoon) and various v- 
platoons. The u-platoons are formed when a bunch of ve- 
hicles escapes from the jam. As argued in Sec. IIII CI the 



first car escapes after an average waiting time T w = j-rji ■ 
Due to anticipation, with probability 1— R the second car 
can move in the same time-step, and so on. Therefore the 
average number of cars escaping in the same time-step is 
given by 



i 

R' 



(8) 



These cars form a w-platoon of length I where the value of 
v depends on the parameter region as discussed above. 
Since the average waiting time for the escape of a car 
is T w , the average distance between two w-platoons is 
Ax c — vT w — i—ji ■ 

To calculate the flow in the congested branch, we again 
neglect the transition regions and assume that only one 
jam with N vehicles and n w-platoons with a total num- 
ber of N v cars are present. Then we have No + N v = N 
with N v — nl. Furthermore Nq + N v + nAx c — L where 
Nq and N v are the total lengths of the platoons and nAx c 
is the total space between the platoons. These relations 
yield 

1 , , N v Ax r , . 

l = -(A + nAx c ) = p+-^--^. (9) 

The average velocity of the vehicles in the congested 
branch is v = Using ©, this implies for the flow 



cong — — 

l-R 



NN V 

v = 1 

L N V 



P) 



AX C 



R 



(10) 



Note that this result is independent of the velocity v of 
the platoons! It is in excellent agreement with the sim- 
ulation data, justifying e.g. the assumption made about 
the transition regions. 



E. Stability regions 

For fixed a we now can estimate the stability region 
pi < p < p2 for the mixed states. At the lower bound- 
ary density p\ the number of cars N v in the u-platoon 
vanishes. From (0 one has p\Axf — 1 = which yields 



Pi 



1 - R 



Vf - v + (1 - R) 



(11) 



The upper bound p 2 is not determined by the condition 
N v = N, i.e. all cars belong to the w-platoon. This would 
correspond to the density p = 1. In fact, the instability 
of the mixed state occurs earlier. At the density 



P2 



R(v + R-2) + l 



(12) 



the flow of the mixed branch becomes larger than that 
of the congested branch, see (|10|l and therefore (at least 
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for random initial conditions) the flow of the congested 
branch is observed. However, our simulations have given 
indications for hysteresis effects and metastability in the 
large density regime. We will discuss these in more de- 
tail in a future publication 27]. For v — the upper 
transition density becomes p2 = 1, independent of R, 
consistent with the observation (Fig. |3J) that the mixed 
region for v — extends up to the maximal density. 

Since p 2 has to be larger than pi, this yields an ad- 
ditional condition for the stability of the branches. It is 
easy to check that pi < p 2 if 

(1-R)v f >v. (13) 

This is just the condition obtained in Sec. IIII CI from 
computer simulations. 

Summarizing, a mixed region with u-platoons can only 
exist for l/(2(v + l)) <a< l/(2v) and R satisfying (|T3)) . 
If these conditions are fulfilled, f-platoons occur in the 
density interval p\ < p < p 2 where p\ and p 2 are given 
by l(TT|l and (fEl)l . respectively. 

IV. SUMMARY AND CONCLUSIONS 

Forecasting the impact of different anticipation 
schemes plays a essential role in real traffic flow in or- 
der to propose automated traffic alternatives. In this 
paper we have introduced and investigated a modifica- 
tion of the NaSch model to better capture reactions of 
the drivers intended to keep safety on the highway. As a 
result, an anticipation parameter a € [0,1], that allows 
to determine a velocity-dependent safe braking distance 
of the precedent vehicle, is included. The addition of this 
parameter proves to be useful to describe different traf- 
fic situations of non-automated, automated, and mixed 
traffic. 

Simulation results presented here for homogeneous 
drivers, corresponding to a cautious estimation of the 
preceding car velocity (large a), are in excellent agree- 
ment with the shape of the empirical fundamental dia- 
gram. 

On the other hand, simulation results for driving 
schemes associated to intermediate-levels of anticipation 
with a from 0.13 to 0.5 (homogeneous drivers), exhibit 
phase separation in a certain density regime into a free- 
flow region and so-called w-platoons. In these dense pla- 
toons vehicles move with the same velocity v and have 
vanishing headway. The velocity v of the platoon is de- 
termined by the level of automation. 

This platoon formation observed in a mixed regime 
plays an important role in Automated Highway Systems 



(AHS) to increase the vehicular capacity. Therefore, re- 
sults obtained help to elucidate the effects of anticipation 
coded in a. The maximum flow and the density interval 
for free-flow regime go with the inverse of a: smaller 
values of a, (greater estimation of the precedent car ve- 
locity) imply larger flows and larger density interval for 
free-flow regime. This is also in accordance with, for ex- 
ample, the use of certain anticipation strategies to exploit 
the knowledge of the velocity of the precedent vehicle and 
so, reducing the distance among vehicles, increasing the 
capacity, and the density interval for free-flow regime. 

Moreover, the analysis of the speed variance of individ- 
ual vehicles also indicates the importance to anticipation 
effects. Results indicate that level of anticipation should 
be determined based on the density. The highest-level of 
anticipation should be considered before the correspond- 
ing maximum density for free-flow regime is reached. Af- 
ter that maximum density, larger density requires lower- 
level of anticipation. Therefore, this selection not only 
produces the highest flow, but also the lowest standard 
deviation, and so, efficiency is the highest. 

The considerations in this paper show the flexibility 
of the CA approach to more complex traffic flow prob- 
lems. A simple and natural modification of the rules of 
the NaSch model to consider different driving schemes 
allows us to describe the formation of coherently mov- 
ing platoons observed in some anticipation schemes. We 
think that the results presented here are relevant to es- 
tablish suitable levels of safety and automation not only 
for AHS, but also in real traffic. We stress that although 
in this paper the model is simulated in a single-lane on a 
ring, it is possible to apply it to complex highway topolo- 
gies in a satisfactory way |27|. 

Apart from its practical relevance for traffic problems, 
our work also shows interesting physical aspects. The 
model suggested here exhibits various kinds of phase sep- 
aration phenomena. At intermediate densities, phase 
separation into a condensed (u-platoon) and a non- 
condensed (free-flow) phase can be observed. In con- 
trast to most other models of driven diffusion, the con- 
densed phase is moving coherently for v > 0. At high 
densities an even more surprising state is found that ex- 
hibits phase separation between different condensates, a 
non-moving one (v = 0) and several coherently moving 
platoons (v > 0). To our knowledge such a behavior 
has not been observed before. It would be interesting 
to study these phases in more detail, especially since re- 
cently some progress in the understanding of phase sep- 
aration in driven diffusive models has been made |28| . 
Work in this direction is currently in progress |27| . 
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